Procesamiento de Datos Municipales

A new article created using the Distill format.

Federico Córdoba https://blackfriedrich.github.io/spatialdataportal/blackfriedrich.html (OFUT-FAU-UNT)https://www.observatoriofau.com.ar/
2024-08-20

Un * al final indica que se ha alojado un dato en el ordenador como producto resultante.

Carga de librerías

Unidad de trabajo

Importación de la data de Municipios y Comunas

Importación:

baseurl_iegwfs <- "http://vport.herrera.unt.edu.ar:8894/geoserver/wfs?"
request_iegwfs <- "request=GetFeature&service=WFS&version=2.0.0&typeName=espacio01:tuc_juris&outputFormat=application/json"

query_iegwfs <- paste0(baseurl_iegwfs,request_iegwfs)
jurisdicciones_df <- st_read(query_iegwfs) #'df' por dataframe
Reading layer `OGRGeoJSON' from data source 
  `http://vport.herrera.unt.edu.ar:8894/geoserver/wfs?request=GetFeature&service=WFS&version=2.0.0&typeName=espacio01:tuc_juris&outputFormat=application/json' 
  using driver `GeoJSON'
Simple feature collection with 121 features and 4 fields
Geometry type: MULTIPOLYGON
Dimension:     XY
Bounding box:  xmin: -66.17672 ymin: -28.02735 xmax: -64.48182 ymax: -26.05595
Geodetic CRS:  WGS 84

Visualización:

mapview(jurisdicciones_df, label = jurisdicciones_df$name, color = "darkgreen", col.regions = "green", alpha.regions = .05)

Intersección de la unidad municipio *

Revisamos la estructura de los datos

str(jurisdicciones_df)
Classes 'sf' and 'data.frame':  121 obs. of  5 variables:
 $ id      : chr  "40" "41" "75" "79" ...
 $ juris   : chr  "SAN ANDRES" "LOS BULACIOS Y LOS VILLAGRA" "MUN. LULES" "SAN FELIPE Y SANTA BARBARA" ...
 $ depto   : chr  "CRUZ ALTA" "CRUZ ALTA" "LULES" "LULES" ...
 $ codigo  : chr  "8154" "" "0077" "8364" ...
 $ geometry:sfc_MULTIPOLYGON of length 121; first list element: List of 1
  ..$ :List of 1
  .. ..$ : num [1:855, 1:2] -65.2 -65.2 -65.2 -65.2 -65.2 ...
  ..- attr(*, "class")= chr [1:3] "XY" "MULTIPOLYGON" "sfg"
 - attr(*, "sf_column")= chr "geometry"
 - attr(*, "agr")= Factor w/ 3 levels "constant","aggregate",..: NA NA NA NA
  ..- attr(*, "names")= chr [1:4] "id" "juris" "depto" "codigo"

Revisamos la lista de nombres de las Jurisdicciones, a partir de la columna “juris”:

sort(unique(jurisdicciones_df$juris))
  [1] ""                                    
  [2] "7 DE ABRIL"                          
  [3] "ACHERAL"                             
  [4] "AGUA DULCE Y LA SOLEDAD"             
  [5] "ALPACHIRI Y EL MOLINO"               
  [6] "ALTO VERDE Y LOS GUCHEA"             
  [7] "AMAICHA DEL VALLE"                   
  [8] "AMBERES"                             
  [9] "ANCAJULI"                            
 [10] "ARCADIA"                             
 [11] "ATAHONA"                             
 [12] "BENJAMIN ARAOZ Y EL TAJAMAR"         
 [13] "BUENA VISTA"                         
 [14] "CAPITAN CACERES"                     
 [15] "CEVIL REDONDO"                       
 [16] "CHOROMORO"                           
 [17] "CIUDACITA"                           
 [18] "COLALAO DEL VALLE"                   
 [19] "COLOMBRES"                           
 [20] "DELFIN GALLO"                        
 [21] "EL BRACHO Y EL CEBILAR"              
 [22] "EL CADILLAL"                         
 [23] "EL CERCADO"                          
 [24] "EL CHAÑAR"                           
 [25] "EL MANANTIAL"                        
 [26] "EL MOJON"                            
 [27] "EL MOLLAR"                           
 [28] "EL NARANJITO"                        
 [29] "EL NARANJO Y EL SUNCHAL"             
 [30] "EL POLEAR"                           
 [31] "EL PUESTITO"                         
 [32] "EL SACRIFICIO"                       
 [33] "EL TIMBO"                            
 [34] "ESCABA"                              
 [35] "ESQUINA Y MANCOPA"                   
 [36] "ESTACION ARAOZ Y TACANAS"            
 [37] "GARMENDIA"                           
 [38] "GASTONA Y BELICHA"                   
 [39] "GDOR. PIEDRABUENA"                   
 [40] "HUASA PAMPA"                         
 [41] "LA ESPERANZA"                        
 [42] "LA FLORIDA Y LUISIANA"               
 [43] "LA RAMADA Y LA CRUZ"                 
 [44] "LA TRINIDAD"                         
 [45] "LAS CEJAS"                           
 [46] "LAS TALAS"                           
 [47] "LEON ROUGES Y SANTA ROSA"            
 [48] "LOS BULACIOS Y LOS VILLAGRA"         
 [49] "LOS GOMEZ"                           
 [50] "LOS NOGALES"                         
 [51] "LOS PEREYRAS"                        
 [52] "LOS PEREZ"                           
 [53] "LOS PUESTOS"                         
 [54] "LOS RALOS"                           
 [55] "LOS SARMIENTO Y LA TIPA"             
 [56] "LOS SOSA"                            
 [57] "MANUEL GARCIA FERNANDEZ"             
 [58] "MANUELA PEDRAZA"                     
 [59] "MEDINA"                              
 [60] "MONTE BELLO"                         
 [61] "MONTEAGUDO"                          
 [62] "MUN.  JUAN BAUTISTA ALBERDI"         
 [63] "MUN.  LA COCHA"                      
 [64] "MUN.  TAFI VIEJO"                    
 [65] "MUN.  TRANCAS"                       
 [66] "MUN.  YERBA BUENA"                   
 [67] "MUN. DE AGUILARES"                   
 [68] "MUN. DE ALDERETES"                   
 [69] "MUN. DE BANDA DEL RIO"               
 [70] "MUN. DE BELLA VISTA"                 
 [71] "MUN. DE BURRUYACU"                   
 [72] "MUN. DE CONCEPCION"                  
 [73] "MUN. DE FAMAILLA"                    
 [74] "MUN. DE GRANEROS"                    
 [75] "MUN. DE LAMADRID"                    
 [76] "MUN. DE MONTEROS"                    
 [77] "MUN. DE SAN MIGUEL DE TUCUMAN"       
 [78] "MUN. DE SIMOCA"                      
 [79] "MUN. LAS TALITAS"                    
 [80] "MUN. LULES"                          
 [81] "MUN. TAFI DEL VALLE"                 
 [82] "PAMPA MAYO"                          
 [83] "QUILMES Y LOS SUELDOS"               
 [84] "RACO"                                
 [85] "RANCHILLOS Y SAN MIGUEL"             
 [86] "RIO CHICO Y NUEVA TRINIDAD"          
 [87] "RIO COLORADO"                        
 [88] "RIO SECO"                            
 [89] "RUMIPUNCO"                           
 [90] "SAN ANDRES"                          
 [91] "SAN FELIPE Y SANTA BARBARA"          
 [92] "SAN IGNACIO"                         
 [93] "SAN JAVIER"                          
 [94] "SAN JOSE DE LA COCHA"                
 [95] "SAN PABLO Y VILLA NOUGUES"           
 [96] "SAN PEDRO DE COLALAO"                
 [97] "SAN PEDRO Y SAN ANTONIO"             
 [98] "SANTA ANA"                           
 [99] "SANTA CRUZ Y LA TUNA"                
[100] "SANTA LUCIA"                         
[101] "SANTA ROSA DE LEALES Y LAGUNA BLANCA"
[102] "SARGENTO MOYA"                       
[103] "SOLDADO MALDONADO"                   
[104] "TACO RALO"                           
[105] "TAPIA"                               
[106] "TENIENTE BERDINA"                    
[107] "VILLA BELGRANO"                      
[108] "VILLA DE CHICLIGASTA"                
[109] "VILLA DE LEALES"                     
[110] "VILLA PADRE MONTI"                   
[111] "VILLA QUINTEROS"                     
[112] "YANIMA"                              
[113] "YERBA BUENA"                         

Seleccionamos la unidad que nos interesa:

municipio_v <- "MUN. DE FAMAILLA" #'v' por value

filtramos

municipio_df <- jurisdicciones_df %>% 
  filter(juris == municipio_v)

visualizamos

mapview(municipio_df, label = municipio_df$name, color = "darkgreen", col.regions = "green", alpha.regions = .05)

Creación de bbox municipio *


Objetos OSM

Importación de objetos OSM (bbox): Calles_bboxmunicipio * AAAA

Importación de objetos OSM (bbox): Amenities_bbox municipio * AAAA


Intersección de objetos OSM (Municipio): Calles_municipio *

Intersección de objetos OSM (Municipio): Amenities_municipio *


Objetos alojados en el ordenador

Importación e intersección de capa parcelaria (bbox): parecelas_bboxmunicipio *

Importación e intersección de capa parcelaria (municipio): parecelas_municipio *


Importación e intersección de capas edificación (bbox): parecelas_bboxmunicipio *

Importación e intersección de capas edificación (municipio): parecelas_municipio *


Visualizaciones:

GGplot: s/bbox | nube de puntos de edificaicón + estructura vial + parcelario + aménitis. Vista Online en blog

GGplot: s/bbox | poligonos de edificación + estructura vial + parcelario + aménitis. Descarga Online o PDF *

Leaflet o Mapviwe: s/Municipio | poligonos de edificación + estructura vial + parcelario + aménitis. Vista Online en blog